UNCLASSIFIED 


AD_  667  642 

METHOD  OF  CALCULATING  TRANSIENTS  IN  GAS  TURBINE 
POWER  PLANTS 

V.  V.  Selin 

Foreign  Technology  Division 
Wright-Patterson  AFB,  Ohio 

14  September  1967 


Processed  for .  .  . 

DEFENSE  DOCUMENTATION  CENTER 
DEFENSE  SUPPLY  AGENCY 


FOR  FEDERAL  SCIENTIFIC  AND  TECHNICAL  INFORMATION 


U.  *.  DEPARTMENT  Of  COMMERCE  /  NATIONAL  BUREAU  OF  STANDARDS  /  INSTITUTE  FOR  APPLIED  TECHNOLOGY 


AD  66  7 


FTD-MT-24 -155-67 


FOREIGN  TECHNOLOGY  DIVISION 

<D 


METHOD  OF  CALCULATING  TRANSIENTS  IN  GAS  TURBINE 
POWER  PLANTS 

by 

V.  V.  Selin 


Distribution  of  this  document 
is  unlimited.  It  may  be 

_  released  to  the  Clearinghouse 

kepfodwCfd  by  the 

CLEARINGHOUSE 
lor  Federal  Sc»ont»fic  &  lccbnitjJ 
Information  Springfield  Va  22151 


Department  of  Commerce,  for 
<sa^l^J;c^Ui€^generj^jub^^^ 


/7 


This  document  is  *  machine  translation  of  Russian 
text  which  has  bean  processed  by  the  AN/3SQ-l6()W-2) 
Machine  Translator,  owned  and  operated  by  the  United 
States  Air  Force.  The  machine  output  has  been  post- 
edited  to  correct  for  major  ambiguities  of  meaning, 
words  missing  from  the  machine's  dictionary,  and  words 
out  of  the  context  of  meaning.  The  sentence  word 
order  has  been  partially  rearranged  for  readability. 
The  content  of  this  translation  dees  not  indicate 
editorial  accuracy,  nor  does  it  indicate  USA?  approval 
or  disapproval  of  the  material  translated . 


FTP-MT-  g^-155-67 


EDITED  MACHINE  TRANSLATION 


METHOD  OF  CALCULATING  TRANSIENTS  IN  GASS  TURBINE 
POWER  PLANTS 

By:  V.  V.  Selin 

English  pages:  12 

SOURCE:  Izvestiya  Vysshikh  Uchebnykh  Zavedeniy. 

Energetika.  (News  of  Institutions  of 
Higher  Learning.  Power  Engineering). 
No.  12,  1963,  pp.  55-6'+. 


UR/0143-6  3-000-012  _ TP7501279 


THIS  TRANSLATION  IS  A  RENDITION  OR  THI  0RWL 

HAL  TOR II ON  TUT  WITHOUT  ANY  ANALYTICAL  OR 
EDITORIAL  COMMENT.  STATEMENTS  OR  THEORIES 

AOVOCATSO  OR  IMPLIED  ARE  THOSE  OR  THE  SOURCE 

AND  DO  NOT  NECESSARILY  REELECT  THE  POSITION 

TRANSLATION  DIVISION 

OR  OPINION  OR  THE  RORIISN  TECHNOLOGY  DL 

ROREION  TECNN0L06Y  DIVISION 

VISION, 

•R.AER,  OHIO. 

FTP- MT-  24-155-67 


Date  14  Sept.  W  67 


ITIS  INDEX  CONTROL  FORM 


01  Acc  Nr 

TP7501P7H 


68  translation  Nr 

FTD-MT- 24 -155-67 


«  x  Ref  Acc  Nr 
AP401440S 


97  Header  Clas|63  Clas  164  Control  Markings 


1  94  Ex pans i 


UNCL 


UN  CL, 


Reel/Frane  Nr 

ixj. _ ojteiL 

lien  40  Ctry  Info 
UR 


02  Ctry 

03  Ref 

04  Yr 

05  Vol 

06  Iss  j 

07  B.  Pg.  45  E.  Pg . 

10  Date 

UP 

loll £ _ 

63 

000 

012 

0055  0064 

NONE 

METODIKA  RASCHETA  PSREKHODNYKH  REZHIMOV  GAZOTURBINNYKH  USTANOVOK 
09  English  Title 

•METHOD  of  calculating  transients  in  gas  turbine  power  plants 

43  Source 

IZVESTIYA  VYSSHIKH  UCHEBNYKH  ZAVEDENIY.  ENERGETIKA (RUSSIAN) 


42  Author  98  Document  Location 


SELIN,  V.  V, 


16  Co-Author 

NONE 


16  Co-Author 
NONE  _ 


16  Co-Author 


47  Subject  Codes 

21 


39  Topic  Tags  1 

gas  turbine,  power  plant,  gas  compressor 


NONE _ _ 

16  Co-Author 

NONE 

abstract:  A  method  of  finite  differences  for  calculating  the  transient 
characteristics  of  a  gas-turbine  power  plant  is  offered;  the  method  is 
claimed  to  be  "simpler  and  more  reliable  as  compared  to  existing  methods." 
The  following  formulas  are  developed:  rate-of-flow  in  the  turbine;  same 
in  the  compressor;  compression  in  the  compressor;  expansion  in  the 
turbine;  balance  of  total  compression  ratio  and  total  expansion  ratio; 
rate-of-change  of  gas  quantity  in  the  capacity  (regeneration  case) ;  heat 
exchange  In  the  regenerator;  power  balance  In  the  compressor  shaft; 
compressor  efficiency;  turbine  efficiency.  The  above  formulas,  some 
additional  turbine  formulas,  and  ship-propulsion  formulas  are  tabulated. 
Peculiarities  in  the  calculation  of  a  load-drop  case  are  discussed.  Orlg. 
art.  has:  ?  figures,  33  formulas,  and  1  table.  English  Translation: 

12  pages. 


Fo** 

FTD  rtt  67  0-90 


U.  S.  BOARD  ON  GEOGRAPHIC  NAMES  TRANSLITERATION  SYSTEM 


Block 

Italic 

Transliteration 

Block 

Italic 

Transliteration 

A 

a 

A 

a 

A, 

a 

P 

p 

P 

P 

R,  r 

B 

6 

£ 

6 

B, 

b 

C 

c 

c 

e 

S,  s 

B 

■ 

B 

I 

V, 

V 

T 

t 

T 

m 

T,  t 

r 

r 

r 

I 

G, 

g 

y 

y 

y 

y 

U,  u 

a 

a 

n 

0 

D, 

d 

4> 

♦ 

* 

P.  r 

E 

e 

E 

i 

Ye 

,  ye;  E,  e* 

X 

X 

X 

» 

Kh,  kh 

)K 

M 

JK 

3tC 

Zh 

,  zh 

u 

u 

u 

H 

Ts,  ts 

3 

a 

3 

t 

z. 

z 

H 

H 

H 

V 

Ch,  ch 

H 

M 

H 

u 

I, 

i 

U1 

tu 

m 

i u 

Sh,  sh 

n 

A 

R 

a 

Y, 

y 

lU 

111 

m 

“f 

Shch,  shch 

K 

K 

K 

K 

K, 

k 

1> 

% 

b 

1 

II 

J1 

n 

n 

A 

L, 

1 

bl 

u 

bi 

y 

Y,  y 

M 

K 

M 

M 

M, 

m 

b 

w 

b 

a 

i 

H 

M 

H 

H 

N, 

n 

9 

a 

3 

a 

E,  e 

0 

0 

0 

0 

o. 

o 

JO 

K) 

!0 

*0 

Yu,  yu 

n 

n 

n 

n 

P, 

P 

a 

n 

ft 

a 

Ya,  ya 

*  ye  initially,  after  vowels,  and  after  t,,  b;  e_  elsewhere. 
wFien  written  as  8  in  Russian,  transliterate  as  yK  or 
The  use  of  diacritical  marks  is  preferred,  but  such  marks 
may  be  omitted  when  expediency  dictates. 


FTD - MT - 24 - 1 r  5 -6  7 


i 


METHOD  OF  CALCULATING  TRANSIENTS  IN  GAS  TURBINE 
POWER  PLANTS 


Engineer  V.  V.  Selin 

Kaliningrad  Mechanical  Institute 
of  Fishing  Industry  and  Economy 


For  a  well-founded  selection  of  a  gas  turbine  power  plant 
[GTUj(rry)  layout,  during  its  designing  of  great  value  is  a  detailed 
calculated  analysis  of  characteristics  in  transitional  modes 
(transients),  for  which  a  sufficiently  simple  and  exact  method  of 
calculation  Is  necessary.  In  this  article  a  method  of  calculation 
is  offered  which  is  more  simple  and  reliable,  as  compared  to  the 
existing  methods. 

Under  variable  conditions,  the  number  of  Independent  parameters 
for  any  layout  of  a  GTU  is  equal  to  the  number  of  combustion  chambers 
[1].  During  transitions  of  a  GTU,  the  number  of  independent  parameters, 
simply  determining  this  or  that  mode,  is  equal  to  the  sum  of  the 
number  of  working  combustion  chambers  and  the  number  of  rotors 
In  the  layout.  As  an  example  let  us  examine  the  two-shaft  layout 
with  a  low  pressure  power  turbine  [TND](THH)  (Fig.  1).  If  for  this 
layout  we  assign,  for  example,  the  temperature  beyond  the  combustion 
chamber  and  the  number  of  rotor  turns,  then  we  will  simply  determine 
the  mode  of  the  power  plant. 

Data  available  in  literature  permit  examining  the  layout  of  a 
high  pressure  air  duct  as  one  which  Is  discrete  [3,  5],  i.e.,  to 
consider  that  the  given  capacity  is  concentrated  in  one  point  —  beyond 
the  regenerator  (Fig.  2). 

One  may  assume  that  establishing  the  gas  parameters  before  the 
turbine  and  after  it  occurs  instantly  [2,  31,  which  is  equivalent 
to  disregarding  the  gas  volumes  in  the  flow-through  part  of  the 
turbine.  Here  the  working  substance  can  be  considered  as  being 
inertlaless.  This  assumption  will  agree  with  the  results  of  experi¬ 
mental  investigations  [3]. 
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Fig.  1.  Diagram  01'  a  two-shaft  ITY 
’with  independent  low  pressure 
power  turbine.  TBJ1  —  high  pressure 
turbine;  TH^  —  low  pressure  turbine; 
K  —  compressor;  PJ  —  regenerator; 

KC  —  combustion  chamber;  FI  —  reduc- 
tor;  1,  2...  -  indices  of  the 
parameters  of  working  substance 
at  a  given  point. 


When  determining  time  t  during  the  transition  from  one  established 
mode  to  another  the  parameters  of  these  conditions  and  character  of 
load  on  power  shaft  are  assumed  to  be  well-known. 

The  problem  is  solved  with  the  help  of  the  method  of  finite 
differences.  One  may  assume  that  the  entire  transient  process 

consists  of  a  number  of  unsteady  modes 
of  the  power  plant,  each  of  which  is 
included  in  a  certain  small,  finite 
interval  of  time  At^. 

As  the  condition  which  disturbs 
the  balanced  state  we  assume  the 
assigned  law  of  change  of  gas  tempera¬ 
ture  beyond  the  combustion  chamber 

T-j  =  f (G  T) ,  taking  into  account  the 

dynamic  characteristics  of  the  control 
system.  Here  and  below  by  a  dash  we 
have  designated  the  relative  magnitudes, 

i.e.,  Gt  =  G-j/G.pQ  —  the  flow  rate  ratio  through  the  turbine  in  a 

given  mode  to  the  flow  rate  in  the  nominal  mode.  Numerical  index 
at  parameters  T  and  p  corresponds  to  the  designations  in  Fig.  1. 

For  determining  the  parameters  of  each  unsteady  mode,  it  is 
necessary  to  solve  a  system  of  equations  which  characterizes  the 
work  of  all  basic  ITy  elements. 

Equation  of  the  flow  rate  for  a  turbine  in  general  has  the  form 

/  (Pjoi  Pm'.  T 2<r  nt><  Pi'.  Pi'.  T 3;  n), 

where  p-j  and  p^  are  the  pressure  before  and  beyond  the  turbine;  is 

the  temperature  of  gas  before  the  turbine;  n  is  the  number  of 
revolutions . 

For  multistage  turbines  it  is  possible  to  use  the  Stodol  equation, 
and  for  low  stage  turbines  when  determining  flow  rate,  we  can  use 
the  method  of  I.  V.  Kotlyar  [1].  Using  the  flow  rate  equation,  for 
the  assigned  law  of  change  of  gas  temperature  before  the  turbine 
Tj  =  f ( G„  ) ,  it  is  possible  to  construct  a  characteristic  of  the 

turbine  p^  =  f(V- 


Fig.  2.  Diagram  of  high 
pressure  gas-air  duct: 

P  —  regenerator;  KC  — 
combustion  chamber;  V  — 
given  capacity  (tank). 
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Equation  of  flow  rate  for  a  compressor  is  assumed  to  be  assigned 
in  the  form  of  a  universal  characteristic. 


Equation  of  compression  in  the  compressor 


H. 


1) 


(1) 


where  is  the  specific  work  of  compression  in  the  compressor; 

Hjj  is  the  isoentroplc  heat  drop  in  the  compressor;  is  the  air 

temperature  at  the  entrance  to  the  compressor;  e  is  the  compression 
ratio;  r)R  is  the  internal  efficiency  of  the  compressor;  k  is  the 

index  of  isoentropy. 


Air  temperature  beyond  the  compressor  is 

r, = r,  -f 

cprm 

Equation  of  expansion  for  a  turbine  is 

”  Wjti,  ss  n,cfT  3  ^  I - 


(2) 

(3) 


where  Lj.  is  the  specific  work  of  expansion  in  the  turbine;  HT  is  the 
isoentroplc  heat  drop  in  the  turbine;  r^T  is  the  internal  efficiency 
of  the  turbine;  cp  is  the  specific  heat  capacity  of  gas;  <p  is  the 
expansion  ratio. 


Qas  temperature  beyond  the  turbine 

r__  t  rT 

*  —  /  j - - . 

cr 


(4) 


and  Ht  are  determined  by  the  thermal  diagram  [1], 

In  certain  cases  one  should  consider  the  nonstationary  heat 
exchange  in  the  flow-through  part  of  the  turbine. 


The  heat  which  is  imparted  by  the  gas  to  the  metal  (by  the 
metal  to  the  gas)  during  the  time  At^, 


SQ 


1  f  it _ r  > 

3600  '  1  e‘  '• 


(5) 


where  a  is  the  average  (along  the  flow-through  part)  coefficient  of 
heat  transfer;  F  is  the  surface,  washed  by  the  gas;  T_  and  T  are, 

i  CT 

correspondingly,  the  average  temperatures  of  gas  and  the  walls  of 
the  blade  channels. 
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A  change  of  gas  temperature 

*Ttr 


_  &Q 
G, f- A-,  ' 


(6) 

(7) 


A  change  of  temperature  of  the  channel  walls 

AT  -  - 
?CT  CM<rM  ' 

where  G„  and  c„  are,  correspondingly,  the  weight  of  metal, 

MM 

participating  in  the  heat  exchange,  and  its  specific  heat  capacity. 

As  a  result  the  temperature  of  gas  beyond  the  turbine  will  be 
determined 


Ti  “  Tj  —  AT,.  (8) 

A  change  of  heat  drop  in  the  turbine  as  a  result  of  the  heat 
exchange  of  gas  with  metal  will  be  proportional  to  the  change  of 
temperature 


(9) 

The  influence  of  thermal  accumulation  in  the  metal  of  the 
combustion  chamber  and  the  gas  conduits  can  be  disregarded  due  to 
the  low  coefficients  of  heat  transfer  and  the  low  metal  content  of 
the  combustion  chamber.  In  many  cases  it  Js  also  possible  not  to 
consider  the  heat  exchange  in  the  flow-through  part. 

Equation  of  balance  of  the  general  compression  ratio  and  expan¬ 
sion  is 

ToCui  T~  » *oCiu’  (10) 

where  £  is  the  loss  factor  of  pressure. 

Equation  of  the  speed  of  change  of  the  quantity  of  gas  in 
capacity^  In  case  of  the  presence  of  regeneration  and  the  gas-air 
capacities  connected  with  it,  flow  rate  through  the  turbine  in 
unsteady  modes  will  be  different  than  the  flow  rate  through  the 
compressor.  This  difference  in  flow  rates  AG^;  at  the  end  of  any 

assigned  Interval  of  time  Ax^  will  be  directly  proportional  to  the 

increase  of  pressure  p  and  Inversely  proportional  to  the  increase 
of  gas  temperature  T  in  the  tank.  As  a  result, 

AG,.,  ■>  (Jlliirji.  _  ,  )  (11) 

V  V-i)  T<  /•  V  1 

where  index  i  pertains  to  the  desired  mode,  while  (i  -  1)  pertains 
to  the  preceding  mode;.  Ry  is  time  of  capacity,  i.e.,  the  time  in 

which  the  tank  V  will  be  emptied  in  the  presence  of  flow  rate  G^,. 


I 

*  I 


•I 


1 


k 


1 


With  a  sufficient  degree  of  accuracy,  the  temperature  of  gas 
in  the  tank  can  be  determined  as  the  temperature  of  air  beyond  the 
regenerator  =  T^.  Value  can  be  found  by  the  trial-and -error 

method,  by  the  characteristics  of  the  compressor  and  the  turbine. 

■j 

Equation  of  heat  exchange  in  the  regenerator.  As  has  been 
already  noted  by  many  authors  [2,  4],  the  accumulation  of  heat  in 
the  metal  of  a  regenerator  renders  a  very  significant  influence  on 
transitions  of  a  TTY.  This  should  be  considered  when  determining 
the  air  temperature  beyond  regenerator  T^. 

The  heat  which  is  imparted  by  the  gas  to  the  regenerator  wall 
(by  the  wall  to  the  gas), 


^  -  •^  (Tr  -Tcr  ). 


(12) 


the  heat,  imparted  by  the  regenerator  wall  to  the  air  (air  to  wall), 

-7>  (13) 

A  change  of  thermal  potential  of  the  wall 


(14) 


In  equations  (12)  and  (13).  and  aB  are  coefficients  of  heat 

transfer  from  the  gas  to  the  wall  and  from  the  wall  to  the  air,  which 
in  the  transient  process  can  be,  with  a  sufficient  degree  of  accuracy, 
taken  as  depending  only  on  the  flow  rates  of  gas  and  air,  if  we  were 
to  consider  the  low  variability  of  coefficients  of  thermal  conduction 
X  and  coefficients  of  viscosity  u: 


*r  ;  :  *r0  Cr ; 


(15) 


where  m  and  n  are  indices  of  degrees  which  enter  into  equation 

Nu  =  a  Rem  and  depend  on  the  type  of  heat  exchanger;  F  is  the  surface 
of  the  heat  exchanger;  Tct,  Tr,  TB  are  averaged  temperatures  of  the 

wall,  gas  and  air. 


Gas  temperature  beyond  the  regenerator  is 


r,  =  7v 


4<?r 


Of  c  p  At, 

Air  temperature  beyond  the  regenerator  is 


r^r.-p 


c«  *>  •»  *, ' 

A  change  of  wall  temperature  of  the  regenerator 


(16) 

(17) 


ifcr  - 


Gu  * 


(18) 


where  gm  and  cm  are»  correspondingly,  the  weight  of  metal, 
participating  in  the  heat  exchange,  and  its  specific  heat  capacity. 

Value  ATqt  is  usually  low,  and  in  many  cases  the  temperature 

of  the  wall  of  the  heat  exchanger  in  trahsient  processes  can  be 
considered  as  being  constant. 

Equati on  of  balance  of  powers  of  the  compressor  shaft.  Excess 
moment  of*  the  compressor  shaft 

AAf,  =  AfTa.i  —  A/*  —  Af,p.  (19) 

where  MTB^,  Mk>  M  are  moments  of  the  turbine,  compressor  and 
friction. 

A  change  of  revolutions  of  the  compressor  shaft 

in,  =  AAf.lv  (20) 

where  I  is  the  moment  of  inertia  of  the  compressor  shaft. 

Equation  of  balance  of  powers  of  the  power  shaft.  Excess 
mome nt  of  the  power  shaft 


AA/e  =  Afina  —  Ma.  —  Mtp.  (  21 ) 

A  change  of  revolutions  of  the  power  shaft 

( 2 2) 

where  Ic  is  the  moment  of  inertia  of  the  power  shaft  (including  the 

reductor  and  the  consumer),  applied  to  revolutions  of  the  power 
turbine.  Mfl3  is  the  moment  on  the  consumer  shaft,  which  can  be 

determined  by  characteristics,  taking  into  account  the  delay  of  the 
control  system.  For  example,  for  a  variable  pitch  screw  ( BPUI)  it 
is  sufficient  to  have  dependences 


R  =  f(v), 


applied  to  revolutions  of  the  power  turbine,  and  a  program  of  screw 
pitch  variation. 

Here  H/D  is  the  pitch  ratio;  v  is  the  speed  of  vessel  in  knots; 
R  is  water  drag  on  vessel  motion;  P_  is  total  thrust  of  the  screws. 
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A  change  of  speed  of  a  vessel  for  the  given  interval  of  time 
At^  will  be  determined  from  equation 

=  (23) 

where  m  is  vessel  mass  with  additional  mass  of'  the  water. 

Internal  efficiency  of  the  compressor  can  be  found  by  its 
characteristic . 

Efficiency  of  the  turbine  can  be  calculated  by  well-known 
dependences  qQJ1  =  f(x)  in  the  variable  mode,  if  x  =  ~  is  defined 
as  the  given  ratio  of  speeds  0 


x  —  xtn  )/  . 


After  determining  all  the  parameters  of  the  i-th  mode,  we  should 
set  as  our  goal  the  time  interval  At^+^  and  find  parameters  of 

the  (i  +  l)-th  mode.  The  calculation  should  be  conducted  up  to 
achievement  of  the  parameters  of  assigned  conditions.  An  example 
of  the  calculation  sequence  is  represented  in  Table  1,  where  two 
unsteady  modes  are  determined  during  reception  of  a  load  for  the 
ITy  (Fig.  I),  working  on  BPUI.  The  law  of  temperature  change  before 
the  turbine  is  accepted  as 


Tt~  const  1  098’  K. 

which  indirectly  considers  the  accumulation  of  heat  in  the 
regenerator. 

According  to  calculated  data  it  is  possible  to  construct 
dependence  n  =  f("r)  and  the  line  of  modes  of  the  compressor  in  the 
transient  process.  Time  of  the  transient  process  t  will  be  determined 
by  the  sum  of  time  intervals  accepted  in  this  calculation 

t  =  2  It,.  (  ) 

i 

The  given  method  can  be  applied  not  only  to  the  simplest,  but 
also  to  any  complicated  layouts  of  TTy.  In  the  presence  of  layouts 
of  ny  of  several  compression  ratios,  division  of  compression  ratios 
by  steps  should  be  made  by  proceeding  from  the  condition  of  their 
joint  work  [1] 


7"““  \/  I—  •  (25) 

where  B  is  the  parameter  of  flow  rate,  equal  to  —1/  —  ;  j  is  the 

*  t  % 

ordinal  number  of  the  compression  ratio. 

Here 

*  =  ....  (26) 
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Table  1 


Designation  of  value 


Interval  of  time 


Number  of  revolutions 
of  compressor  shaft 


Relative  number  of  revolu¬ 
tions  of  compressor  shaft 


Compression  ratio 


Formula 


t --  we  assign 


«—  we  assign 


11 


mrnen- 1  Value 
sion 


ijCw< 


Relative  flow  rate 
through  the  compressor 


Relative  flow  rate 
through  the  turbine 


Pressure  beyond  the  KHJJ 


Compression  ratio  (check) 


Pressure  beyond  the  TB.Q 


low-  «) 

(by  characteristic 


Pi^  Pf  f(G,) 


p» 


Expansion  ratio  of  TBJ1 

Pi 

?•  --7 

P) 

- 

1.51 

1.58 

Heat  drop  of  TB.Q 

11  x  **  /  (?i: 

)  ketlAx 

29.5 

33.2 

Coefficient 

X  |  *3  rtn 

| 

0.312 

.0.32 

Efficiency  of  TRfl 

5i  -H*x) 

Power  of  TB.Q 

N,  -  4,(9  G,GtW,»„ 

kV 

1176 

1525 

Moment  of  TBH 

Af,  •=  973 

nu 

m 

430 

505 

Efficiency  of  compressor 


Heat  drop  of  compressor 


Air  temperature  beyond 
the  compressor 


Power  of  the  compressor 


Moment  of  the  compressor 


Excess  moment  of  the 
compressor  shaft 


v-by  characteristic 


H*  =-  /(•; 


A’k-  4.19  0,0,^- 


y\ih  =-  m,  -  uK  -  w, 


Table  1  (continued) 


Designation  of  value 


Number  of  revolutions  of 
power  shaft 


Relative  number  of  revolu¬ 
tions  of  power  shaft 


l 

Expansion  ratio  of  TH,Q 

l 


'emperature  before  THJ1 


Heat  drop  of  TH# 


Coefficient 


Efficiency  of  TH.H 


Power  of  THU 


Moment  of  TH.H 


Speed  of  the  vessel 


Relative  screw  pitch 


Moment  of  the  screw 


Excess  moment  of 
power  shaft 


Change  of  revolutions 
of  power  shaft 


Thrust  of  the  screw 


Water  drag  to 
vessel  travel 


vessel  speed 


Temperature  beyond  the  THJ[ 


Fo rmula 


•*"«  ^  it7„  i" 


'*<  —  nC(l_ I) 


Dimen- 


Table  1  (continued) 


* 

Formula 

Dimen¬ 

sion 

Value  | 

Designation  of  value 

t 

t+l 

Temperature  beyond  the 
regenerator1 

bbbehh 

fl 

807 

790 

Specific  gravity  of  gas 
in  the  tank 

•(  —  34 1 ,5 

k«/»3 

0.608 

0.77 

Time  of  capacity  (tank) 

1In  this  case  tempera 

plified  formula,  since  it 
and  Rv;  r  is  the  degree  of 

n  __Ri£L 

Ku  —  r 

ture  Tc  is  determined 

5 

is  used  only  when  dete 
regeneration. 

•  |  3,07 

by  a  sim- 

rmining  7 

3 

In  case  of  a  IT/  layout  without  regeneration,  the  calculation 
is  simplified.  For  layouts  with  intermediate  heat  feed,  it  Is 
necessary  to  assign  a  law  of  change  of  gas  temperature  in  the  Inter¬ 
mediate  combustion  chamber. 

Peculiarities  of  calculating  fTY  layouts  for  load-drop.  In 
order  to  determine  the  dynamic  stoking  of  ny  revolutions  during 
full  load-drop  (for  example,  cutoff  of  a  generator  from  the  network), 
it  Is  necessary  to  know  the  value  of  time  lag  of  the  control  system 
Tp  (the  time  which  has  passed  from  load-drop  up  to  cessation  of  fuel 
feed).  After  cessation  of  fuel  feed,  pressure  of  gas  before  the 
turbine  will  be  determined  by  the  pressure  of  gas  in  the  tank.  For 
the  case  of  isoentropic  expansion,  the  connection  between  the 
parameters  of  gas  in  the  tank  will  be  expressed  as 


PI  Polo*, 

where  7  is  the  specific  gravity  of  gas. 

By  differentiating  this  equation  we  will  obtain 


To 


dr. 


In  the  presence  of  small  changes  of  p,  ratio 


I. 


Thus , 


dv 


■  dp  — . 


FTD  -  MT  -  - 1 >  -6  7 
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A  change  of  weight  quantity  of  gas  in  the  tank  can  be  determined 
as 

da  -  Vd;. 
or 

da  —  (C,  —  G,)  th. 

If  we  consider  that  in  the  normal  mode  GTq  =  GK0,  then,  crossing 

to  finite  Increments,  after  certain  transformations  we  can  determine 
the  relative  change  of  pressure  in  the  tank  as  a  result  of  expansion 
for  the  time  interval  Ax^. 


Ap 


kJ±< 


(27) 


One  can  assume  that  after  cessation  of  fuel  feed  into  the 
combustion  chamber,  gas  temperature  before  turbine  at  any  moment 

of  time  is  equal  to  the  exit  temperature  from  the  tank.  In  turn, 
gas  temperature  in  the  tank  at  the  end  of  time  interval  Ax^  will  be 

determined  by  the  temperature  of  displacement  TCJ|  of  gas  which  is 

anew  proceeding  to  the  tank  from  the  regenerator,  during  the  time 
Ax^  with  the  remaining  gas  in  the  tank  and  the  change  of  gas 

temperature  as  a  result  of  expansion  AT 


Tt~Ten-lT. 


(29) 


If  during  interval  of  time  Ax1  into  the  tank  quantity  of  gas 

GkAXi  with  temperature  proceeded  from  the  regenerator,  and  from 

the  tank  quantity  of  gas  (^Ax^^  with  temperature  departed  into 

the  turbine,  then  from  the  equation  of  heat  balance,  if  one  were 
to  accept  the  process  of  mixing  as  being  instantaneous,  we  can 
obtain  the  temperature  of  mixing  in  the  tank 


j" _ (ff,  *"  7*3  ■* 


A  change  of  gas  temperature  in  the  tank,  as  a  result  of 
isoentropic  expansion,  will  be  expressed  by  equation 

or  _j_ 

■=.  k—  I  ~  •  , 

dT  =  — -  p  dp. 

Considering  that  in  the  presence  of  small  changes  of  pressure  value 
~p  *  as  I,  by  crossing  to  finite  increments,  we  can  obtain 


•ip. 


(?0) 


PTD-MT- 24 -159-67' 
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For  any  Intermediate  mode  the  capacity  time  is 


(31) 


According  to  the  data  of  detailed^alculation  for  load-drop, 
it  is  possible  to  construct  dependence  n  =  f(x)  and  to  determine 
the  value  of  dynamic  stoking  of  turns  AETmax* 

It  is  possible  to  prove  that  for  two  PTy,  carried  out  according 
to  the  same  layout,  with  identical  gas  parameters  but  having 
different  power  and  different  design  execution  (different  moments 
of  inertia  of  the  rotors  and  number  of  revolutions)  this  equality 
is  true 


} 


* 


~~rf  =  17  —  const.  (32) 

lnl 

where  n  is  the  parameter  of  pickup. 

For  two  similar  PTY  of  different  power,  dynamic  stoking  of 
revolutions  will  be  identical 

^*inm  =  IflJm***  (33  ) 

Equations  (32)  and  (33)  permit  (in  necessary  cases)  producing 
a  reduction  of  dynamic  PTY  characteristics. 
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